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Abstract 



The Feynman-Kac path integration problem was studied in the worst case setting by 
Plaskota et al. (J. Comp. Phys. 164 (2000) 335) for the univariate case and by Kwas 
and Li (J. Comp. 19 (2003) 730) for the multivariate case with d space variables. 
In this paper we consider the multivariate Feynman-Kac path integration problem 
in the randomized and quantum settings. For smooth multivariate functions, it was 
proven in Kwas and Li (2003) that the classical worst case complexity suffers from 
the curse of dimensionality in d. We show that in both the randomized and quan- 
tum settings the curse of dimensionality is vanquished, i.e., the number of function 
evaluations and/or quantum queries required to compute an ^-approximation has a 
bound independent of d and depending polynomially on e . The exponents of these 
polynomials are at most 2 in the randomized setting and at most 1 in the quantum 
setting. Hence we have exponential speedup over the classical worst case setting and 
quadratic speedup of the quantum setting over the randomized setting. However, 
both the randomized and quantum algorithms presented here still require extensive 
precomputing, similar to the algorithms of Plaskota et al. (2000) and Kwas and Li 
(2003). 
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1 Introduction 



In this paper we study the multivariate Feynman-Kac path integration prob- 
lem. Multivariate Feynman-Kac path integrals are path integrals over the 
space of continuous functions from M.+ to M. d , equipped with a Wiener mea- 
sure. The multivariate Feynman-Kac path integral is the solution of the initial 
value problem for the multivariate heat equation, see Section 2. This type of 
path integral plays a significant role in many fields, e.g., quantum physics and 
chemistry, differential equations, and financial mathematics. A brief survey of 
research concerning path integrals is contained in the introduction of [15]. 

In this paper we continue the research initiated in [15], where a new algorithm 
for computing Feynman-Kac path integrals was proposed. That paper dealt 
with the univariate case (i.e., with one space variable) and the algorithm 
presented there was based on /^-approximation. An extension of the approach 
of [15] to the multivariate case (with many space variables) was presented in 
[7]. It turns out that algorithms based on L 2 -approximation are no longer 
applicable in the multivariate case. The multivariate case can be solved by 
using uniform approximation as a basic building block. In both papers [7,15] 
the problem was studied in the worst case setting for input functions belonging 
to a class F. Assuming that the uniform approximation problem for the class 
F has worst case complexity of order e~ Q ( F ) for some positive a(¥) it was 
proved in [7] that the number of function evaluations required to compute an 
^-approximation is roughly of the same order e~ a< - ¥ \ We stress that typically 
a(F) depends on the number d of space variables and tends to infinity with d, 
in which case we have the curse of dimensionality. 

In this paper, we consider the multivariate Feynman-Kac path integration 
problem in the randomized and quantum settings. We present algorithms that 
compute an e- approximation and analyze their cost. These algorithms are also 
based on uniform approximation. However, the power of randomization and 
quantum computation permits the improvement of the worst case complexity 
bound 0{e~ a ^). Namely, the number of function evaluations required by the 
randomized algorithm is roughly of order £ - 2a ( v )/( a ( ¥ )+ 2 ) ^ whereas the number 
of function evaluations and queries required by the quantum algorithm is 
roughly of order e _a ( F )/( a ( F ) +1 ) ) see Section 7.1. We stress that the exponent 
of e _1 in the randomized setting is at most 2, and in the quantum setting is 
at most 1. 

In addition to providing the algorithms in the randomized and quantum set- 
tings we also study the complexity of multivariate Feynman-Kac path integra- 
tion in the randomized and quantum settings. As in [7,15], the complexity is 
bounded from below by the complexity of multivariate weighted integration. 
The upper bounds are provided by the costs of the algorithms presented in 
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this paper. 



For the class F of r times continuously differentiable functions we have a(F) = 
d/r, and so the worst case setting suffers form the curse of dimensionality. In 
the randomized setting, the complexity is roughly of order 0(e _2 ^ 1+2r//rf - ) ), 
whereas in the quantum setting it is roughly of order 0(e~ 1 ^ 1+r ^ d ^). In both 
cases the curse of dimensionality is vanquished. We thus have exponential 
speedup over the worst case setting. For d ^> r, we have quadratic speedup of 
quantum complexity over randomized complexity. 



2 Multivariate Feynman-Kac path integration 



The multivariate Feynman-Kac formula is the solution of the initial value 
problem for the heat (diffusion) equation 

dz 

— (u,t) = lAz(u,t) + V(u)z(u,t) for (u,t) E intR d x [0, oo), (1) 
*(u,0)=v(u). (2) 

Here v , V : R d — > R are the initial value function and the potential function, 
respectively. As usual, A denotes the Laplacian. 

The solution z of (1) and (2) is given by the Feynman-Kac formula 

z(u, t) = J v(x.(t) + u) exp (jT V(x(s) + u) ds^j w(dx). (3) 

Here, C is the set of continuous functions x : R + — > M. d such that x(0) = 0. The 
path integral (3) is with respect to the <i-dimensional Wiener measure w, see 
[5,16]. Obviously, (3) holds for functions v and V for which the path integral 
exists. In what follows, we assume that the functions v and V belong to a 
class F for which (3) exists. This class is defined in Section 4. 



3 Problem 



For a given fixed point (u*, t*) Gl^x [0, oo) and arbitrary functions v, V from 
the class F, we want to compute an ^-approximation of the exact solution 

Z v y(u*,t*)0i (3). 

The ^-approximation a v y(u*,t*) is computed by an algorithm A n that uses n 
function values of v and V, i.e., 

a vy (u*, t*) = A n (u*, **, u(ui), . . . , u(u fc ), V(u fc+ i), . . . , V(u n )) . 
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3.1 Worst case setting 

In the worst case setting the error of the algorithm A n is defined as 
e WOI (A n )= sup \z v y(u*,t*)-a vy (u*,t*)\. 

v,V€F 

We want to determine the minimal number 

n wor (e, F) = min{n : 3 A n such that e WOT (A n ) < e} 

of function values that are needed to compute an e-approximation in the worst 
case setting. This setting was analyzed in [7]. 

3.2 Randomized setting 

In this setting we use randomized algorithms and replace the worst case error 
assurance by an expected assurance. A randomized algorithm A n depends on 
a random element u> chosen from some probability space Q. More precisely, 
we compute 

a v y(u*, t*; u) = A n ^ (u*, t*, u^.i), . . . , v(u Ujk ), V(u U!>k+1 ), . . . , Vfu^J) , 

(4) 

with n = E aJ (n w ). This means that we allow a random choice of a mapping A n>ul 
and sample points u^, as well as the number of sample points, whose 
expected value is fixed and equal to n. 

We measure the randomized error of the algorithm A n with respect to the L 2 
norm, i.e., 

e rand (AJ:= sup (EMv,v(u*,t*) - a vy (u*,t*;Lu)) 2 Y /2 . 

v,vgf v y 

As before, we want to determine the minimal expected number of function 
values 

n rand (5, F) = min{n : 3 A n such that e rand (A n ) < e} 
needed to compute an ^-approximation in the randomized setting. 

3.3 Quantum setting 

In the quantum setting we use quantum algorithms with (deterministic or 
randomized) quantum queries and assume that we can also perform function 
evaluations and arithmetic operations on a classical computer. These classical 
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operations are used to prepare an input for a quantum algorithm and to trans- 
form the outcome of a quantum algorithm to an approximation of the exact 
solution. We will be interested in minimizing the total number of quantum 
queries and function evaluations needed to compute an ^-approximation. 

In this section, we give a brief overview of a simplified quantum model of 
computation for continuous problems and describe deterministic and random- 
ized quantum queries. We refer the reader to [1,2,3,8,14,21] for more detailed 
information. 

We first outline a general framework of the quantum setting. Assume that for 
a given class H of input functions / : D — > C we want to approximate the 
solution operator 

S : H -> G, 

with G being a normed space whose norm is denoted by || • We will ap- 
proximate S(f) by a quantum algorithm defined below. 

First, we transform a given input function / G H by using a classical algo- 
rithm P s with s classical function evaluations and obtain 

f = P s (f):D^C. 

Then we use the transformed function / as an input to a quantum algorithm. 

Quantum algorithms are defined as follows. Let C 2 be the two dimensional 
complex vector space. Let Hk = C 2 <8> • ■ ■ <8> C 2 be the fc-fold tensor product of 
C 2 , having dimension 2 fc . Let U n (f) : H n — > H n be a unitary operator of the 
form 

U n (f) — QnQ fQn-1 • • " Ql QfQo, 

with unitary operators Qo, ■ ■ ■ ,Q n an d a quantum query Qj, for some / G 
P S (H). 

The deterministic quantum query Qj is defined as in [2]. Let Tik = 'Hm < 8>T~ik-m 
for some m < k. Then Qj : Ti m ® H-h-m ~^ ® 7~Lk-m is a unitary operator 
of the form 

Qf\x)\y) = \x)\y®(3(f(T(x)))), (5) 

with 

t : {0, . . . , 2 fc - 1} -> D, (3 : C -> {0, . . . , 2 k ~ m - 1} 

and © denoting the addition modulo 2 fc ~ m , see again [2] for a more detailed 
discussion. 

The randomized quantum query is defined in [21]. In this case, Qj = Qj iUJ 
depends on a random element uj and Qf u has the form (5) with r = t w 
depending on uj. This permits the computation of approximate values of / at 
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randomized points. Hence in this case the unitary operator U n depends on a 
random element 00 and has the form 

U n ,Uf) = Qn^ Qf,w Qn„-i ■ ■ ■ QiQf,u,Qo with n = Eriu. 
As usual we assume that the initial state is |0) and we compute 

= U n (f)\0) = Q n QfQn-i ■ ■ ■ QiQfQo\0) 
for deterministic quantum queries and 

=U n „(f)\0) =Q n Q f -, u Qn-i---QiQ f -, u Qo\0) 

for randomized quantum queries. Then we measure the final state and obtain 
an outcome j G {0, . . . , 2 k — 1} with probability 

pf(j) = mj}\ 2 or p ftU u) = m, u \j)\ 2 . 

Knowing the outcome j we compute the final result on a classical computer, 
and the quantum algorithm A n yields 

MfJ) = 00') or A n ^(f,j) = <p w (j). 
for some <fi or <f> u . 

In this paper we will be using quantum algorithms with randomized quantum 
queries. The error of such an algorithm A n is defined as 

e^(A n , P„ S) = sup flE^Eq - A n ^(P s (f), j)\\ 2 G ) 1/2 , (6) 

where E w is the expectation over the probability space Q, and E q is the ex- 
pectation with respect to distribution of the quantum algorithm outcomes. 

Similarly to the other settings, we want to determine the minimal number of 
random quantum queries and classical function evaluations 

n quant( £j jjj = min | s + n . ^ Pg ^ A n SUch that e qUant (A n , P s ) < s} 

needed to guarantee that the error does not exceed e. 

Remark 1 We now briefly comment on the quantum error setting defined 
by (6). Let us concentrate for a moment on the randomness introduced by 
a quantum algorithm, leaving aside randomized queries. So far, the literature 
dealing with continuous problems in the quantum setting has mainly considered 
probabilistic error. That is, instead of taking an expectation with respect to all 
possible outcomes of a quantum algorithm fos E q in (6)), we want an error 
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estimate such that 

(EjS(f)-A n , Ps(fl J 2 G ) 1/2 <e 

holds with a certain (high) probability, for any f G H . Obviously these two 
ways of measuring the error of a quantum algorithm are related. We choose to 
study the average error for simplicity. Moreover, the average error is probably 
more natural when we consider randomized queries. 

The multivariate Feynman-Kac path integration problem in the quantum set- 
ting is defined by taking / = (v, V) with H = F x F and S(f) = z v y(u*, t*). 



4 The function class F 



To assure the existence of the path integral (3), we need to choose a proper 
class of input functions F, see [7]. We assume that 

F = {(f 1 ,f 2 )E¥x¥: ||/i||f<A, ||/ 2 ||f < M (7) 

is a ball of a linear space F x F for some positive /3 2 - 

We make the following assumptions about the linear space F. 

(1) We assume that for every u G M d , the functional L u : F — > R defined by 
L u f = /(u) is continuous, and for arbitrary a,t G R+ we have 

J c ||-^x(t)||Fexp [a ||L X ( S )|| F ^ w{dx) < oo. (8) 

By the Fernique theorem, see e.g., [6], condition (8) holds if there exists 
a < 2 such that ||£ x ||f = 0(||x||°) for ||x|| approaching infinity, see [15] 

for details. Here and elsewhere, ||x|| = \JYh=\ x 1 is the Euclidean norm 
in R d . 

(2) We assume that F is continuously embedded into L OQ (R d ). That is, F C 
L 00 (M, d ) and there exists a positive K such that 

ll/IU^) < K\\f\\ r V/GF. (9) 

This assumption permits us to relate the multivariate Feynman-Kac path 
integration problem to uniform approximation in the worst case setting, 
see again [7]. 

By uniform approximation we mean the approximation of the embed- 
ding operator S : F — > L 0O (M d ), Sf = f in the norm of L 00 (M d ). Let 
n APp( £ > ^) denote the minimal number of function values needed to com- 
pute an ^-approximation in the worst case setting. As we shall see in 
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Section 6.1, uniform approximation also plays a significant role in the 
randomized and quantum settings. 
(3) We assume that 

nl™ P (e,W) = 0{e' a{¥) ) as e -> 0, (fO) 

for some positive number a(¥). 

The linear space F is characterized by the exponent a(¥). Usually a(F) de- 
pends on the smoothness and the number of variables of functions in F, see 
Section 9. 



5 Feynman-Kac formula as a series of multivariate integrals 



In this section we briefly recall some results from [7] which are needed for our 
analysis. 

Without loss of generality we can assume u = in (3). Then we can express 
the path integral as a series of multivariate integrals 

oo 

S(v,V):=z(0,t) = Y,S k+1 (v,V), (11) 

k=0 

where 

k 

S k+1 {v,V)= v(z k+1 )T\V(z i )g k+1 (z 1 ,...,z k+1 )dz 1 ...dz k+1 , (12) 

1=1 

with 

9k+i( z i, ■ ■ ■ , z fe+i) = / fk+i (ti, ■ ■ ■ , t k , t, z 1; . . . , z k+1 ) dti . . . dt k 

J0<ti<---<t k <t 

(13) 

and 

/ fc+ i t k , t, Zi, . . . , z fc+ i) = ((27r) fc+1 *i(* 2 - ti) • • • (t - t k )) ~ d ' 2 

/ 1 / ||zi|| 2 ||z 2 - zi|| 2 ||z fe+ i - z fc || 2 \\ 

xexp V 2 V *i *2-*i "' // ' 

Note that the integral (12) depends on the input functions v and V only 
through the product 



h k+ i{z u . . . , z fc+ i) = u(z fc+ i) ]J ^( z 



i=l 
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and the weight functions g k +i can be computed in advance, albeit with diffi- 
culty. Let us recall also that 



t k 

ls , fc+i|L 1 (]R(fc+i)d) = £j for k > 0. (14) 



6 Approximation of one term of the series 



In this section we present algorithms approximating one term of the series (11). 
To make the notation more clear we define a weighted integration operator 

4+1(7) = / 1W /(zi, • • • , z fc+ i) ^ fc+ i(zi, . . . , z fc+ i) dz 1 ... dz k+1 

where / : R( fc+1 ) d — > R is an integrable function. We can then rewrite one term 
of the series (11) as 

S k+1 (v,V) = I k+1 (h k+1 ). 

In both the randomized and quantum settings, we shall use deterministic uni- 
form approximation of the function h k+ \. To utilize the power of randomization 
and/or quantum computation, we will apply the known technique of variance 
reduction. 



6.1 Variance reduction 



Smolyak's algorithm is a powerful tool for computing an ^-approximation of 

fc+i 
, *• . 

tensor product problems. For h k +i G F ® • • • ® F, Smolyak's algorithm is of 
the form 

n(e,k+l) 

,k+l(hk+l) — ^2 ^fc+l(ti,e,l> • • • ) ti,e,fc+l)Ci,e,fc+l) (15) 
i=l 

for some t i:£ j E M. d and Ci,e,fe+i e L oc (]R( A:+1 ) a! ). It is proven in [7, Lemma 2] 
that 

\\h k+ i - C4,fc+i(/ifc+i)|U 00 (R(*+i)««) < £ \\ v \\v \\V\\w, ( 16 ) 

where 

n( £ ^ + l)< Co ^ 1 + C2 ^J e-^ ¥ \ (17) 

for some q G K. Here a + denotes max{a,0}, the right hand side of (17) is 
defined to be c e~ a ^ when fc = 0. 
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The idea underlying variance reduction idea is as follows. First we compute 

hk+l,e = U £) k+l{hk+l) 

using n(e, k + 1) function values. Then we compute 

Ji(e,fc+1) 

^fc+l(^fc+l,e) — ^fe+l(ti, £ ,l, • • • , ti,£,k+l)h+l(Ci,e,k+l)- 

i=l 

Observe that the functions Ci,e,fc+i do not depend on the input functions v and 
V so the integrals ifc+i(Ci,e,fc+i) can be precomputed. 

We stress that hk+i, £ and I k +i{hk+i,e) are deterministic. We will use random- 
ized or quantum algorithm to approximate the multivariate integrals 

Ik+l{hk+l — hk+l, £ )- 

Since the error depends on the norm \\hk+i — hk+i t£ \\L oc (K( k + 1 ') d )> which is small, 
we can do this efficiently. We present the details in the following two sections. 

6.2 Randomized algorithm 

To make formulas simpler we define 

fk+l,e = hk+i — hk+l, £ - 

We use the randomized algorithm of the form 

C™Vl>> V) = h + l{h k+1 , £ ) + Q^ +1 Jfk + l,e). (18) 

Here 

osSn^(/) = -EM>) ( 19 ) 

m 3=1 

denotes the classical Monte Carlo algorithm with m randomized sample points. 

Randomized sample points are chosen with respect to the density 
9k+i/ \\gk+i\\ L 1 (R( k + 1 ) d ) which is indicated by the random parameter uj e VL. 

Using the well known error formula for the classical Monte Carlo algorithm, 
we conclude that 

(E,(/ w fei)-ci + i,i^^)) 2 ) 1/2 

( ( - rand - \ 2 V /2 

= (I k+ l(fk+l,e) - Q r ™k+l,u>(fk+l,e) ) J 

= _L (Var(/ fc+li£ )) 1/2 , (20) 
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with 

Var(/ fe+li£ ) = h+i{fk+i,e) ~ (4+i (jW)) • 
Clearly, from (16) and then from (7), (14) we get 

(Var(/ fe+li£ )j < - ||/ fe+ i, £ || Loo(K(fc+ i) d) < £ £j < £ fcj ■ 

(21) 

This yields the error estimate 

(E„(W^ + i) - C£i>> ^)) 2 ) 1/2 < -4 ^ff^ ( 22 ) 



and the total number of function evaluations is 

n(e,k + l)+m. (23) 

6.3 Quantum algorithm 



The structure of our quantum algorithm is similar to randomized one, having 
the form 

Cw,>' V ) = h+l(hk+l,e) + QmTk+l,dMk+he), ( 2 4) 

with, as before, fk+i,e = hk+i — hk+i,e- Here, we use a quantum algorithm 
QmTk+ii w hh /t randomized quantum queries, that approximates the classical 
Monte Carlo algorithm (19). In [4] the problem of approximating 

-E/(^>) 

m i=i 

was analyzed for Boolean functions /. Using the technique of reducing the 
summation problem for bounded real functions to the summation problem for 
Boolean functions as in [2], we see that a result similar to that of [4] holds. 
From [4] and (16) we conclude that 

x 2\ V2 

I TO \ " ' 

E q (-E/(^)-«XU.(/^, 

= o (^-ll/fc+i, £ ||L 00 (R( fc + i ) d )) = o (^r^ 1 ^) • 

By integrating over Q, we obtain 

2 V /2 ^ ( £ A#2* fc> 



E w E q 



Q T m*k+l,w(fk+l,e) - QZ^k+l^{fk+l,e) J - O \- ^-^y j . (25) 
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The total number of queries and function evaluations is 

n(e, k + 1) + k, 

We stress that this number does not depend on m, which is only used for the 
definition of the Monte Carlo algorithm. 

We now estimate the total error as 

1/2 



< ( E^ (/fc+l(/fc+l,e) - Qm^k+l^ifk+lr 



1/2 



This, by (22) and (25), yields 

(^C«(V,) - C£, t+1 ,>,iO) 2 ) 1/2 = o (-£= + £ . 

Letting m = k 2 we get the error bound 

(E.E q (/ fc+1 (/^ + i) - CSU + i>« ^)) 2 ) 1/2 = O (26) 

using 

n(e,k + l) + K (27) 
function values and quantum queries. For the sake of convenience we denote 

I quant /quant i 2 

<P e ,K,k+l*, = <P e ,m,K,k+lfi, Wlth m = K . 



7 Complete algorithms 



Based on the previous two sections we are ready to present algorithms com- 
puting an e-approximation of multivariate Feynman-Kac path integral. We 
approximate consecutive terms of the series 



S(v,V) = J2Sk + i(v,V) 

k=0 



by the algorithms 



irand J.quant 
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with the accuracies e™+f and in the corresponding settings being 

rand _ 2/(q(F)+2) ^' quant _ l/(q(F)+l) ^ /an\ 

k+1 (3^t k 2 k +^ k+1 (3 1 (3$t k 2 k + 2 { ] 

and the number of randomized sample points rrik+i and quantum queries Kk+i 
being 

(29) 



m k+1 



£ -2a(F)/(a(F)+2) 



£ -a(F)/(a(F)+l) 



The final forms of randomized and quantum algorithms approximating S , (v, V) 
are 



and 



fc=0 

^y-quant 

$ quant ( y) = £ ^quant 



where the finite integers i\T5 and and A^ uant will be determined in the next 
section. 



7.1 Error analysis 

From (22), (26) and (28), (29), it is easy to check that the error bounds 

(EUh + i(h k+1 ) - ^, mfc+1>fc+1 >, V)f) < (30) 

and 

1/2 



(E„E q (I k+1 (h k+1 ) - 0^ +i>fc+1 > V)f) < JL. (31) 



hold. 



It is also easy to see that we need to approximate only a few terms. Indeed, 
for k approaching infinity, we have ££+i anc ^ ^fc+T* a ^ so tending to infinity. 
Note that by (9), (22) and (26) we see that for 

,-rand quant 

^±L > K k+1 and > K k+ \ 



with K being the embedding constant in (9), the deterministic zero algorithms 
provide sufficient accuracy. Thus, we need to use the algorithms 4> r ^ d mfc+i fc+1 u 

and , , only for k = OQnfe- 1 )). Hence, we get iV5 and = 0(\n(e^)) 

and iV; quant = C^lnfc- 1 )). 
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The bounds (30) and (31) yield 

(E„(S(v,V) - ^(v,V)) 2 ) 1/2 < e, (32) 
(E w E q (S(t;, V) - $^ a >, V)) 2 ) 1 ' 2 < e. (33) 

This means that the algorithms $^ and and $^ uant compute ^-approximations of 
the multivariate Feynman-Kac path integral in the randomized and quantum 
settings respectively. 



7.2 Number of function values and quantum queries 



In this section we derive estimates on the numbers of function values and 
quantum queries n($™ nd ) and n($^ uant ) of the algorithms $^ and and <J>^ uant . 
By the bounds (17), (23) and (27) we get the obvious estimates 

mi/£ fe+1 



n ($r d ) = O ( ( + E ( ci + ^-^f^j 



k=i 



x 



' a ftk-i-knk+l \ a ( ¥ ) \ 



and 

, 1 / quant \ (a(F)+l)fc 
mi / £: fc+l \ 



n 



x ^l^ 2fe+1 ^ ( ) + ^g-afFJ/CaCFJ+l) ] _ 

We can now use an argument similar to that in the proof of [15, Theorem 1] 
to show that 

oo / i •. /_rand\ (.a(F)+l)k / fl ^fc.fcofc+lX a ( F ) 



fc=i 



oo 



Inl/erfV am+1) * flhlW2 k+1 ' 



a(F) 



for all 5 > 0. Thus we finally get 

n ($- nd ) = O ( e -2«CF)/(a(F)+2)-«) , (34) 

n ($ £ quant ) = O ( £ -««/(a(F)+i)-^ (35) 

for all 5 > 0. 
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8 Complexity of multivariate Feynman-Kac path integration in 
randomized and quantum settings 

An analysis of the complexity of the multivariate Feynman-Kac path integra- 
tion in randomized and quantum settings is quite similar to the one presented 
in [7] and [15]. We only point out essential differences. 

8.1 Lower bounds 

Lower bounds for our problem complexities are provided by the complexities 
of multivariate weighted integration problem. By this problem we mean an 
approximation of the integration operator / : F — > R define by 

1(f) = (2ntT d/2 f /(u) exp(-||u||/(2t*))du V / e F. 

J» d 



Consider a randomized algorithm A™ that uses n function values and ap- 
proximates the integration operator I . We say that this algorithm computes 
an e-approximation of the weighted integral if 

(E w (/(/)-A^ d ) 2 ) 1/2 <e V/eF. 

We denote by rag^e, F) the minimal number of function values needed to 
compute an e-approximation in the randomized setting. 

Consider a quantum algorithm A quant that uses n randomized quantum queries 
and approximates the operator /. We say that A^ &nt computes an ^-approximation 
of the weighted integral if 

(E w E q (/(/) - ^ u r) 2 ) V2 < e V/ e F. (36) 

We define n < Y^ t (s, F) as the minimal number of quantum queries needed to 
compute an e-approximation. 

As in [7], we can reduce multivariate Feynman-Kac path integration to mul- 
tivariate integration with a Gaussian weight by taking 1^ = 0, since S(v, 0) = 
I{v). Moreover, (7) and (9) imply that 

n\^(e,¥)<n^ d (e,¥), 
n^\e,¥) < n quant (e,F). 
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8.2 Upper bounds 



Obvious estimates on the complexity of the multivariate Feynman-Kac path 
integration are provided by the cost of the algorithms derived in Section 7. 
Thus, by (34) and (35) we get 



n rand (e,F) = O ( e -2a(F)/(«(F)+2M) , 
^quant^j^ = q ^- Q (F)/(a(F)+l)-^ 



for all 5 > 0, where a(F) is the exponent of the uniform approximation problem 
complexity for the space F containing the class F, i.e., 



nX°p r P (5,F)=0(e^ F )), 



see also Section 4. 

From the previous two sections we can see when the randomized and quantum 
algorithms proposed in this paper are almost optimal. This is the case for the 
classes of input functions for which randomized and quantum complexities of 
the integration problem defined in Section 8.1 are of orders e~ 2a ( F )/( a ( F )+ 2 ) and 
£ -a(F)/(a(F)+i) respectively. 



9 Examples 



In this section we present two examples of function classes F satisfying the 
assumptions from Section 4 and compute lower and upper bounds of of the 
complexities of the multivariate Feynman-Kac path integration. 



Weighted Sobolev space 

We use one of the results from [19], which relates the complexity of the ap- 
proximation of functions defined over a finite domain to the complexity of the 
weighted approximation of functions over the whole space M. d . Let 

F=j/:r*^R: f e C r (R d ) and ||/|| f := ]T ||/ (a) || Loo(Rd) < oo 

I 0<|a|<r 

For simplicity, we consider a weight function p : M. d — > R + , given by 

p(z) = exp(-||z|| 2 ) VzGl d 
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which decays exponentially. By [19] there exists an algorithm 



that computes a weighted ^-approximation of the function / G F, i.e., 

IK/-W)p||L.(R-)<e||/llF. 

with 

n = (e- d/r ) . 

Let 

¥ = {f:R d ^R: f/peW and ||/|| F := ||//p|| f < oo}. 

We can use the algorithm U§ to construct an algorithm U £ approximating 
functions from F. Indeed, define 

U J = Y,f(he) a ie, where a it 



i=i 



Then for f p = f / p, we have 

||/ - U £ f\\ L , Rd) = ||(/, - U>f p ) p\\ Lac(Rd) < e ||/J f = e 



as claimed. We have to check the three remaining conditions which are to 
be satisfied by F, namely, the continuity of function evaluation as well as 
conditions. It is easy to see that for / G F and z 6 l d we have 

/(z) < \\f\\ Locm < Wf/Ph^) < II/IIf 

and so function evaluation is continuous. Conditions (8) and (9) follow imme- 
diately from this continuity. 

The algorithms $^ and and $^ uant compute an ^-approximation of the multivari- 
ate Feynman-Kac path integration problem for the class F with the number 
of function evaluations and/or quantum queries roughly 0(s^ 2 ^ 1+2r ^ d ^) and 
0(e~ 1 ^ 1+r ^), respectively. However, the factors appearing in the big O no- 
tation depend on d and this dependence is exponential, see Sections 7.1 and 
7.2. For d ^> r the exponents Idjid + 2r) for the randomized algorithm and 
d/(d + r) for the quantum algorithm are close to 2 and I. In fact, the orders 
2 and 1 can be obtained by the use of the classical Monte Carlo algorithm 
(without variance reduction). Then, the factors multiplying e~ 2 and e~ l are 
independent of d for the class F so the curse of dimensionality present in the 
worst case setting (see [7]) is indeed broken when we switch to the randomized 
or quantum settings. 
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To obtain lower bounds on the multivariate Feynman-Kac path integration 
problem in the class F we may switch to the integration problem as in Sec- 
tion 8.1. We observe that this integration problem is not easier than the uni- 
form integration over the unit cube 

1(f) = f /(z) dz 

by taking functions with support [0, l] d . It is known that the uniform integra- 
tion problem has the randomized complexity Q(e~ 2 ^ 1+2r ^ d ^), see [12,13], and 
quantum complexity 0(£ -1 /( 1+r / d )), see [21]. This shows that the algorithms 
(jjrand an j ^quant are r^g^y optimal for the class F . 



Periodic functions 

This example was considered in [7]. We repeat all details for the reader's 
convenience. Following [17] we consider the class F of 27r-periodic functions 
/ : [0, 2n] d — > R d satisfying the condition 

V/GF Vj = l,...,d 3^GL 0O ([-27T,27r] (i ) 

1 r 2w 

/(x) = — / (pj(x u ...,Xj-t,...,x d ) F r (t) dt, (37) 

Z7T JO 

where r > and 

00 / r 7r\ 

F r {t) = l + 2j2 k ~ r cos { kt ~~2j- 

The norm in the class F is defined as 

1 d 

ll/llf = j IIV 5 illL 00 ([-2 7 r,27r]d), 

where the ipj are functions from the representation (37) of the function /. In 
[17], there is a linear algorithm U £ that computes a uniform ^-approximation 
of functions from the class F, i.e., that 

\\f-U £ f\\ L ^ [QMd) <e\\f\\ w V/GF, 

with the cost of order e~ d l r . 

Denote by F the class of functions / : M. d — > R that are periodic extensions of 
functions from F. Let ||/||f : = ||/|[o,27r] d llF- Obviously, problem of the uniform 
approximation for the class F can be obviously solved using the algorithm 
mentioned above with the same cost as for the class F. Similarly to the previous 
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example, we have to check the three conditions of Section 4. It is easy to see 
that for / G F, z £ M. d , and arbitrary j G {1,2, ... ,d} we have 

/(z) < ||/||l°°(]R<») = ll/l[0,27r] d ||L°°([0,27r] d ) < C \\ (fj || L=°([-27r,27r] d ) 5 

with C = (27T)- 1 J 27r |F r (t)| dt. Hence 

|/(z)|< ll/llx-^j <C ||/|| F 

and so function evaluation is continuous. The remaining conditions follow 
immediately. 

Thus the algorithms $^ and and $q uant ; based on the algorithm U £ described 
above, compute an ^-approximation of the multivariate Feynman-Kac path in- 
tegral with a number of function evaluations and/or quantum queries roughly 
0(E~ 2 /( 1+2r / d }) and 0(e~ 1 '( 1+r / c v), respectively. Using an argument similar to 
that of the previous example, we conclude that the algorithms $^ and and $^ uant 
are roughly optimal for the class F. 
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